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Asstract. In this paper we discuss the algorithmic generation of Penrose- 
Stairs. Furthermore we will work out several proves, type-classifications and 
investigate numerical and geometrical behaviour of those staircases under dif- 
ferent preconditions. Graphics and coding-examples in different programming- 
languages are being added to illustrate our subjects. 


1. INTRODUCTION 


In the year 1937 the swedish graphic artist Oscar Reutersvard created the draw- 
ing of a fours-sided staircase, whose stairsteps lead all the way down, ending finally 
exactly at their starting position again, to form a 2D-parallel projected closed ge- 
ometrical object that looks like being constructable in 3D-space, but factual is an 
illusion that’s impossible to exist in 3D-space. Later in the year 1958 two eng- 
lish mathematicians, Lionel Penrose and his son Sir Roger Penrose independently 
discovered and made popular the “impossible staircase”. 


FIGURE 1. The original Penrose-Staircase. 


Though many interesting articles has been written since and lots of beautiful 
graphics are created still nowadays, there seems to be the lack of a serious investi- 
gation how to mathematically and algorithmically obtain a perfect looking Penrose- 
Staircase in parallel projection graphics. What we are looking for in addition, is 
the set of all existing Penrose-Staircases S' = {51, 52, 53,... Sn}. 


Date: 2022. 
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2. COORDINATE SYSTEM 


The well-known original image of the Impossible-Staircase (1) is actually iso- 
metric, although back in 1958 the picture had been squeezed in y-size (perhaps to 
obscure its isometric character). Later we will see and discuss its exact derivation. 
In order to obtain a perfect parallel-projected Impossible-Staircase, we start with 
a basic isometric coordinate-system. We determine, that in flat isometry the coor- 
dinates are made out of equilateral triangles and we declare the sidelength of one 
of those triangles being one single unit u=1. 


FIGURE 2. Isometric grid with triangle sidelength of one single unit. 


Since the original Penrose-Staircase doesn’t fit perfectly in integer isometric- 
coordinates and since we don’t have an algorithm yet, we use a slightly different 
staircase that we easily found by trial-and-error: 


FIGURE 3. Impossible staircase that fits perfectly into an isometric grid 


We assert the count of stairsteps for each side a=6 b=2 c=2 d=6, and the side- 
length of each single stairstep 1=3 isometric units, following written as a/b/c/d(1) 
e.g. 6/2/2/6(3). We define a,b,c,d € Nyo. Also it’s important that every side of 
the staircase a,b,c and d contains at least 2 stairs to assure, that going up and down 
from every single stair leads into two different directions. We define: co > a < 2, 
co>b<2,~ >c¢K<2,~>d<2. 
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3. WALKING DOWN THE STEPS 


Next we need a walking-routine and an ISO-position-cursor, so we define the 
walkable directions we can go from an ISO-point to another like this: 


U 


UL UR 


1/2 


DL DR 


D 


FIGURE 4. Walkable ISO-directions. 


Note that height h in an equilateral triangle is determined by cos(30°)=0,866025403... 
To project our walking functions into cartesian coordinates, which can be repre- 
sented by a computer-screen, we define them as follows: 


U(n):y=ytn 


D(n):y=y-n 


1 
UR(n):c=a+nh; y=yt+n-; 


2 
1 
DR(n):a2=a4+nh; y=y ns: 
1 
DL(n) : a =a — nh; y=y— nz; 


1 
UL(n): «2 =a — nh; y=ytns; 
Using that functions walking downstairs, we get our path which starts at point 
S, going all the way down and around the staircase until it finally ends at point S 


again, defining the shape of our 6/2/2/6(3) staircase: 
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FIGURE 5. The walking path from startpoint S down the staircase 
and back to S. 


Firstly, by extracting and joining the walking-formulas for x we get x = ah + 
bh —ch— dh which is x = h(a+b—c—d). Our x/y-path must start and end at (0,0) 
so we can say 0 = h(a+b—c—d). Because of distance h > 0 always, that equation 
only gets zero if (a + b — c— d) is zero, so our final equation for the x-direction is: 


0=a+b—c-d 


Secondly, by extracting and joining the walking-formulas for y we get 


7 (=) a—b—c-d+4 


Same as above with x, we need to set y to zero, in order to end up at the starting- 
point S again, and then we solve for the single-stairside-length 1, so we get a more 
expressive value in result than just zero: 


a+b+c+d-4 


(= 


Division by zero occurs here for staircases with all sides having the same count 
of steps: a=b=c=d, or if a =cAb=d. To avoid this, we declare a 4 c. 

Also, our result, the length 1, has to be always greater than zero to make sense, 
thus we can assume that if | <= 0 the staircase is not valid, as well as for | = co. 
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Because of 0=a+b—d—c we can saya+b=c+d. The sum g of all stairs 
around the staircase is 
g=a+b+c+d-4 


therefore we get a+ b= até andc+d= at With this knowledge it’s possible to 
solve for a,b,c and d: 


_g9t4 


ar b 
p= 22", 
+4 
ont -d 
d= i", 


Remembering our formula for lenth 1: 


p etbterd=A 


we see that the numerator is our stairsum g and the denominator is the determinant 


(r): 


With our knowledge so far we can write a computer-program that finds the length 
1, stairsum g and a,b,c,d from 2 up to 100 by brute-force method in FreeBASIC: 
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> written and (c) 2022 by F. Lehr 
print "Searching Penrose-Stairs by Brute-Force-Method" 


dim a as integer 

dim b as integer 

dim c as integer 

dim d as integer 

dim 1 as double ’sidelength of a single stair 

dim x as double 

dim n as integer = 100 ’ max. stairs of one side a,b,c,d. 


Open "Penrose_abcd_2_to_100.csv" For Output As #1 
print #1, "a,b,c,d,numstairs,len," 


for a=2 to n 
for b=2 to n 
for c=2 to n 
for d=2 to n 
x = at+tb-c-d 
if x=0 and a<>c then 
l=(a + b+c+d- 4)/(a/2 - b/2 - c/2 + d/2) 
if 1>0 and 1<10 then 
print #1, ""; a; ","3; bs ","3 c3 ","5 ds "4"5 
Catbtctd-4); ","; 1; "," 
end if 
end if 
next d 
next c 
next b 
next a 


Close #1 
print "press key ..." 
sleep 


We obtain a CSV-Table that we can sort ascending by the stairsum g. Because 
we have chosen a maxlength up to 100 for every side a,b,c,d we can be sure, that 
we will find the complete set of the first 35 Penrose-Staircases. 
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n/| a/b] cl] dj g len 
1} 3] 2] 2] 3] 6 6 
2/3} 3] 2] 4] 8 8 
3] 4} 2] 2] 4] 8 4 
4} 4] 2] 3] 3] 8 8 
5} 3} 4] 2} 5] 10 10 
6/ 4} 3] 2} 5] 10 5 
7} 4} 3] 3} 4] 10 10 
8} 5} 2) 2] 5] 10 | 33333333333 
9} 5} 2] 3} 4] 10 5 

10); 5} 2] 4] 3] 10 10 
11) 3] 5] 2] 6} 12 12 
12); 4} 4] 2] 6} 12 6 
13) 4} 4] 3] 5] 12 12 
14) 5} 3] 2] 6} 12 4 
15) 5} 3] 3] 5] 12 6 
16) 5] 3] 4] 4] 12 12 
17) 6} 2] 2] 6} 12 3 
18) 6} 2] 3] 5} 12 4 
19); 6| 2] 4] 4} 12 6 
20} 6) 2} 5] 3] 12 12 
21} 3] 6] 2} 7) 14 14 
22) 4| 5] 2} 7) 14 7 
23} 4| 5] 3] 6) 14 14 
24} 5| 4] 2) 7| 14 | 4,6666666667 
25} 5| 4] 3} 6) 14 7 
26} 5| 4] 4} 5) 14 14 
27) 6] 3] 2} 7) 14 3,5 
28| 6] 3] 3] 6| 14 | 4,6666666667 
29} 6] 3] 4] 5) 14 4 
30} 6] 3] 5] 4) 14 14 
31} 7] 2} 2) 7) 14 2,8 
32} 7/ 2] 3} 6) 14 3,5 
33 | 7| 2] 4] 5 | 14 | 4,6666666667 
34] 7|/ 2] 5) 4) 14 7 
35} 7] 2] 6} 3) 14 14 
As we can see in the table, the length | is sometimes equal to the stairsum g. 
This can happen only if the determinant r in 1 = £ becomes 1. Also we see, that 


the stairsum g gets exactly a-* times equal to length 1. Noticable is, that length 
1 can be a rational number sometimes! This means, that our ISO-Cursor drifts 
off the grid during it’s walk, but ends exactly where it has started. Hey, we just 
found the original Penrose-Staircase from 1958! It’s the 28th Staircase and has a 
singlestair-sidelength of 4.666666... which is actually ae 


From looking at the table we can easily guess the algorithms for the development 
of a and c as n increments by 1 (see code-listing examples below, function AINC(p) 
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and CINC(p)): 


Aseries = {3, 4, 4, 5, 5, 5, 6, 6, 6, 6, fiaue } 
Cseries = {2, 2, 3, 2, 3, 4, 2, 3, 4, 5, see } 


Now we have enough information to solve our equations above for b = at4 - 


a 


and d = a+ b—c, and length | (remembering our formula for 1) and to write a 
program that calculates the n-th Penrose-Staircase. 


The n-th stairsum is given by: 2n + 4. Thus, if we calculate the amount of Stair- 
cases that share the same n-th stairsum, then the number k of occurring Staircases 
up to the n-th stairsum (inclusive) is: 


b= So flan +4) 


i=l 


4. VARIOUS 


The smallest penrose-staircase is: 81 = 3/2/2/3(6). 


The amount of staircases Q within the n-th stairsum x seems to be: 


1 1 
Q(n) = st ene 


5. CODE EXAMPLES 


Penrose-Staircase-Database CSV-File Generator for PYTHON: 


HHFHHHHEHHHAEHHEHEHHHEEH RARER RHEE ERR HERR SS HH 


# He HH H H FH 


PENROSE -STAIRCASE-GENERATOR v1.0 
This program calculates ratios for 
the "Impossible-Staircase". 

(c) 2022 by F. Lehr 
https://www.ferdinandlehr.de 
ferdinand@ferdinandlehr.de 


# HH HH HH H HF 


HHFHHHHHHHHAEHHEHEHHHRRHEH RRR HEHEHE RHEE S SHH 


import csv 


#Calculates the number of staircases for a given (always even) 


stairsum. 
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#def P(g): 
# return int ((1/8)*g**2-(3/4)*g+1) 


#Calculates the number of staircases for a given (always even) stairsum. 
def P2(g): 

c=0 

p=l1 

k=2 

i=6 

while (i<=g): 

#print("i,g,c",i,g,c) 


c=0 

for j in range(p): 
ct+=1 

#print (c) 

pt=k 

k=k+1 

i+=2 


return c 


#Set of all PStairs to an upper limit : nth Stairsum 
def tz(n): 
r=0 
for i in range(1, n+1): 
yr += P2(2*i+4) 


return int(r) 


#Calculates the nth stairsum 
def NG(n): 
return (n*2+4) 


#Calculates the stairsum g for the nth Pstair 
def SNP(n): 
k=1 
while (tz(k)<n): 
k=k+1 
return NG(k) 


#calculates how many stairsums occur up to a given g (including stairsum for g) 
def SBG(g): 
n=1 
while (NG(n)<g): 
n=n+1 
return n 


#Calcs 1 of the nth PStair. 
def LONP(n): 


10 FERDINAND LEHR 


g=SNP (n) 

#print ("Stairsum =",g) 

rpos=P2(g) 

#print ("Number of stairs with that stairsum p =", rpos) 
m=int((g-4)/2) #calculate how many times d=l for a given g 
#print ("Position within that stairsum range rpos =", rpos) 
#print("m =", m) 

tpos = tz(SBG(g))#n+rpos-1 

#print ("Cursor position total tpos = ", tpos) 

u=m 


for j in range(mt1): 
d=1 #divisor 
1=-1 #length 
for i in range(u): 
l=g/d 
d=d+1 
tpos-=1 
#print ("1 =", 1, "tpos=", tpos) 
if (tpos+1)==n: 
return 1 


u=u-1 


# returns 3,4,4,5,5,5,6,6,6,6,... 
def AINC(p): 
n=1000000000 
c=0 
for i in range(1,n+1): 
for k in range(1,i+1): 
c+=1 
if(c==p): 
return i+2 


#calcs a of the nth pstair 
def A_of_NP(n): 
g=SNP (n) 
rpos=P2(g) 
tpos = tz(SBG(g)-1)+1 
#print ("rpos,tpos,",rpos,tpos) 
return (AINC(n-tpost1)) 


#calcs b of the nth pstair 
def B_of_NP(n): 
a=A_of_NP(n) 
g=SNP (n) 
return int (((g+4)/2)-a) 


# returns 2,2,3,2,3,4,2,3,4,5,2.. 
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def CINC(p): 
n=100000000 
c=0 
for i in range(1,nt+1): 
for k in range(1,i+1): 
#print (i,k) 
ct+=1 
if(c==p): 
return k+1 


#calcs c of the nth pstair 
def C_of_NP(n): 
g=SNP (n) 
rpos=P2(g) 
tpos = tz(SBG(g)-1)+1 
#print ("tpos=",tpos) 
return (CINC(n-tpost1)) 


#calcs d of the nth pstair 
def D_of_NP(n): 
a=A_of_NP (n) 
b=B_of_NP (n) 
c=C_of_NP (n) 
return at+b-c 


#calculate and print the nth pstair 
def PStair_nth(n): 

a=A_of_NP(n) 

b=B_of_NP (n) 

c=C_of_NP (n) 

d=at+b-c 

g=SNP (n) 

1=LONP (n) 

print(a,b,c,d,g,1) 


#calculate and print ist pstair to nth pstair 
def PStairs_to_n(n): 
for i in range(1,n+1): 
PStair_nth(i) 


#calculate and write ist pstair to nth pstair to CSV-file 
def PStairs_to_n_CSV(n): 
f = open(’pstairs.csv’, ’w’, encoding=’UTF8’, newline=’’) 
writer = csv.writer(f,delimiter=’,’) 
for i in range(1,n+1): 
a=A_of_NP (i) 
b=B_of_NP (i) 
c=C_of_NP (i) 
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d=at+b-c 
g=SNP (i) 
1=LONP (i) 
#print(a,b,c,d,g,1) 
#data = [a,b,c,d,g,str(1).replace(’.’, ",")] 
data = [a,b,c,d,g,1] 
writer.writerow(data) 
f.close() 


#Main-Program 

print ("PENROSE-STAIRCASE GENERATOR v1.0") 

n=10 

print ("Calculating the first", n, "Penrose-Stairs") 
print("a b cd sum len") 

PStairs_to_n(n) 
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Viewer for “Context Free Art”: 


ala aa // 
// THE INFINITY-STAIRCASE VIEWER (CFDG) // 
// (c) and written 2022 by F. Lehr // 
// For more info visit https://www.ferdinandlehr.de // 
// github.com // 
// https://www.contextfreeart.org // 
// This code is Creative Commons licensed: // 
// Creative Commons Attribution -NonCommercial - // 
// ShareAlike 4.0 International: // 
// https://creativecommons.org/licenses/by-nc-sa/4.0/ // 
[ [ann eo See ee eo eo eee see // 
CF:: Impure=1 

//----------------- PLAYAROUND-AREA ------------------- // 


COLOR1 = [b 0.007 sat 0.5 h 0] 
COLOR2 = [b 0.5 sat 0.5 h 150] 
COLOR3 = [b 1 sat 0.2 h 200] 


// R = 3,2,2,3,6 // Smallest Staircase 
R = 6,3,3,6,4.66666 // Penrose-Staircase 1958 
// R = 24,6,14,16,5.6 // For Reutersvard and Escher 
// R = 15,2,5,12,3 // Castle of the Wizards 
// R = 17,17,10,24,9.14285 // Large Antique Arena 
// R = 100,66,67,99,9.93939 // Chinese Wall 
// For more ratios visit https://www.ferdinandlehr.de 
/[/-------------------------- 2-2-2 --- 2-2-2 ---- 2-2 ------ // 
U= 1 
H = 0.866025404 
L = R[4] 
A=R [0] 
B=R[1] 
C=R[2] 
D=R[3] 
WH = Dx*2 


startshape PAGE 
shape PAGE { 
PSTAIR[skew 2 6 r -30] 
} 
shape PSTAIR { 
loop i=A-1,-1,-1 [] f{ 


S(L)Ex ((U*0.5*L)*i+(U/2)*i) y (((H*L)-H)*i)] 
} 
loop m=1,D-1,1 [] ¢f 
xx1 = ((U*0.5*L)*(A-1)+(U/2)*(A-1)+(L¥1+U/2)*(B-1)) 


yyl = (((H*L)-H)*(A-1) -(B-1) *H) 
xxxl = (xx1-L*U*0.5*(C-1)+(U/2)*(C-1)) 
yyyi = (yy1-(C-1)*H*(L+1)) 


S(L)Cx (xxx1-L*U*m+(U/2)*m) y (yyy1-H*m) z (-m+A+B+C+D+1)] 
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} 
loop j=1,B,1 [] { 
SCL) Ex ((U*0.5*L)*(A-1)+(U/2)*(A-1)+(L*1+U/2)*j) y (CCH*L) -H)*(A-1)-j*H)] 
} 
loop k=1,C,1 [] { 
xxl = ((U*0.5*L)*(A-1)+(U/2)*(A-1)+(L*1+U/2)*(B-1)) 
yyl = (((H*L)-H)*(A-1)-(B-1)*H) 
S(L) Cx (xx1-L*U*0.5*k+(U/2)*k) y (Cyy1l-k*H*(L+1)) z (A+B+C+D+2)] 
} 
inner_wall2[] 
front_wall[z (A+B+C+D+3) ] 
mid_wall2[] 
right_wall[] 
loop n=2,B,1 [] { 
px = ((U*0.5*L)*(A-1)+(U/2)*(A-1)) + L*U*n + (n-1)*(U/2) 
py = ((CH*L)-H)*(A-1)) -H*(n-1) 
stair_rect_C[x px y py] 
} 
loop o=1,C-1,1 [] ¢f 
px = ((U*0.5*L)*(A-1)+(U/2)*(A-1)) + L¥U*B + (B-1)*(U/2) - L¥U 
py = (((H*L)-H)*(A-1)) -H*(B-1) 
stair_rect_D[x (px-L*U*0.5*0+(U/2)*o0) y (py-(L+1)*H*0)] 


} 
path stair_rect_D { 
MOVETO (0,0) 
LINETO (U*L,0) 
LINEREL (U/2, -H) 
LINEREL (-U*L, 0) 
CLOSEPOLY () 
FILL() [trans COLOR1] 
STROKE(0.1, CF::RoundJoin) [] 
} 
path stair_rect_c { 
MOVETO (0,0) 
LINETO (U*L*0.5,H*L) 
LINEREL (U/2, -H) 
LINEREL (-U*L*0.5, -H*L) 
CLOSEPOLY () 
FILL() [trans COLOR2] 
STROKE(0.1, CF::RoundJoin) [] 
as 
path stairlinesc { 
MOVETO (0,0) 
MOVEREL (L*U*0.5+U*2* A+L+(L/3)+U/2,H*L+H*2* (A-1) -H) 
loop i=0,(B-2),1 [] ¢f 
MOVEREL (U/2, -H) 
MOVEREL (-L*U*0.5, -H*3) 


} 
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MOVEREL (-U/2,H) 
MOVEREL (L*U*0.5,H*L) 
MOVEREL (U*0.5+L, -H) 

} 

MOVEREL (-L-U*1.5, -H*L) 

MOVEREL (U/2, -H) 

MOVEREL (-U*0.5*L, -H*L) 

loop j=0,(C-2),1 [] f 
LINEREL (U*0.5, -H) 
MOVEREL (-L*U*0.5, -H*L) 

} 

FILL()[b 1] 

STROKE(0.1, CF::RoundJoin) [] 


path stairlinesB { 


} 


MOVETO (0,0) 

MOVEREL (L*U*0.5+U*2* A+L+(L/3)+U/2,H*L+H*2* (A-1) -H) 

loop i=0,(B-2),1 [] f 
LINEREL (U/2, -H) 
MOVEREL (-L*U*0.5, -H*3) 
LINEREL (-U/2,H) 
MOVEREL (L*U*0.5,H*L) 
MOVEREL (U*0.5+L, -H) 

} 

FILL()[b 1] 

STROKE(0.1, CF::RoundJoin) [] 


path mid_wall2 { 


} 


MOVEREL (L*U+L*U*0.5+U/2,H*L-H) 
loop n=1,A,1 [] ¢{ 
LINEREL (L*U*0.5,H*L) 
LINEREL (U*0.5, -H) 
} 
LINEREL (-L*0.5*U, -H*L) 
LINEREL (U*WH*0.5, -H*WH) 
LINEREL (-U*0.5*L*(C-2) , -H*L*(C-2)) 
CLOSEPOLY () 
FILL() [trans COLOR2] 
STROKE(0.1, CF::RoundJoin) [] 


path inner_wall2 { 


MOVEREL (L*U+L*U*0.5+U/2,H*L-H) 
loop n=1,A,1 [] f{ 
MOVEREL (L*U*0.5,H*L) 
MOVEREL (U*0.5, -H) 
} 
MOVEREL (-L*0.5*U, -H*L) 
loop m=1,B-1,1 [] { 


15 
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FERDINAND LEHR 


LINEREL (L*U ,0) 

LINEREL (U*0.5, -H) 
} 
LINEREL (L,0) 
LINEREL (U*0.5, -H) 
LINEREL (-L,0) 
LINEREL (U*WH*0.5-(B*U*0.5) , -H*WH+H*B) 
LINEREL (-L*U*(B-2)+U*0.5, -H) 
CLOSEPOLY () 
FILL() [trans COLOR1] 
STROKE(0.1, CF::RoundJoin) [] 


path right_wall { 


} 


MOVETO (0,0) 

loop i=1,(D),1 (] { 
MOVEREL (L,0) 
MOVEREL (-U*0.5,H) 

} 

MOVEREL (L,0) 

loop i=1,C,1 [] { 
LINEREL (L*U*0.5,H*L) 
LINEREL (-U*0.5,H) 

} 

LINEREL (L*U*0.5,H*L) 

loop j=1,cC,1 [] { 
LINEREL (U*0.5, -H) 

} 

LINEREL (U*WH*0.5, -H*WH) 

LINEREL (-U*0.5*L*C, -H*L*C) 

CLOSEPOLY () 

FILL() [trans COLOR2] 

STROKE (CF: :RoundJoin) [] 


path front_wall f{ 


} 


MOVETO (0,0) 

loop i=1,(D),1 [] { 
LINEREL (L,0) 
LINEREL (-U*0.5,H) 

} 

LINEREL (L,0) 

LINEREL (U*WH*0.5, -H*WH) 

LINEREL (-U*D*L,0) 

CLOSEPOLY () 

FILL() [trans COLOR1] 

STROKE(CF::RoundJoin) [] 


path S(1) { 


MOVETO (0,0) 
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LINETO(U*0.5*1,H*1) 
LINETO(U*0.5*1+U*1,H¥*1) 
LINETO(U*1,0) 
CLOSEPOLY () 

FILL() [trans COLOR3] 
STROKE (CF: :RoundJoin) [] 


